Role of delay in the mechanism of cluster formation 
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We study the role of delay in phase synchronization and phenomena responsible for cluster for- 
mation in delayed coupled maps on various networks. Using numerical simulations, we demonstrate 
that the presence of delay may change the mechanism of unit to unit interaction. At weak cou- 
pling values, same parity delays are associated with the same phenomenon of cluster formation and 
exhibit similar dynamical evolution. Intermediate coupling values yield rich delay-induced driven 
cluster patterns. A Lyapunov function analysis sheds light on the robustness of the driven clusters 
observed for delayed bipartite networks. Our results reveal that delay may lead to a completely 
different relation, between dynamical and structural clusters, than observed for the undelayed case. 
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Studying the impact of network topology on dynamical 
processes is of fundamental importance for understanding 
the functioning of many real world complex networks [l[ . 
The dynamical behavior of a system depends on the col- 
lective behavior of its individual units. One of the most 
fascinating emergent behavior of interacting chaotic units 
is the observation of synchronization 0] ■ In general, syn- 
chronization may lead to more complicated patterns in- 
cluding clusters [3|-|5[ . The interplay between underlying 
network structure and dynamical clusters has been the 
prime area of focus for the last two decades @ . Further- 
more, communication delay naturally arises in extended 
systems 0] • A delay gives rise to many new phenomena in 
dynamical systems such as oscillation death, stabilizing 
periodic orbits, enhancement or suppression of synchro- 
nization, chimera state, etc 

In this paper, we study the impact of delay on the phe- 
nomenon of phase synchronized clusters in coupled map 
networks. We investigate the formation of clusters on 
various networks namely, 1-d lattice, small-world, ran- 
dom, scale-free and bipartite networks [TBj . and provide 
a Lyapunov function analysis for bipartite networks to 
explain possible reasons behind the role of a delay on 
synchronized clusters. So far, studies on delayed cou- 
pled dynamical systems mostly concentrated on a global 
synchronized state, except a few recent studies which 
have focused on pattern formation or clustered states 
[1, [H, [H, [I?}- These studies have revealed that delay 
emulates qualitative changes in clustered state, whereas 
mechanism of delayed unit to unit interactions has not 
been investigated so far. 

Previous studies on undelayed coupled systems have 
identified two different phenomena for synchronization 
namely, the driven (D) and the self-organized (SO) Q. 



SO (D) synchronization refer to the state when clusters 
are formed because of intra-cluster (inter-cluster) cou- 
plings. Here, we report that a delay can play a crucial role 
in the formation of clusters as well as the phenomenon 
behind it. The formation of delay-induced synchronized 
clusters may be because of inter-cluster couplings, in- 
stead of coupling between synchronized units [5|, 1 1 61 ] . In- 
troduction of a delay may result in a transition from D 
to SO synchronization or vice versa. Furthermore, our 
studies demonstrate a delay-induced emergence of dy- 
namical phase synchronized D patterns. These patterns 
are stable with time and are dynamical with respect to a 
change in r. A delayed bipartite network leads to a tran- 
sition from SO to D synchronization in an intermediate 
coupling range irrespective of r. 

Here wc take networks with a less average degree 
(Nc ~ N), leading to phase synchronized clusters in- 
stead of a complete synchronized state which usually 
spans all the nodes. We consider a network of N nodes 
and N c connections between the nodes. Let each node 
of the network be assigned a dynamical variable x % , i = 
1,2, ... ,N. The dynamical evolution is defined by the 
well known coupled maps (l8j . 
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Here A is the adjacency matrix with elements Aij tak- 
ing values 1 and depending upon whether there is a 
connection between i and j or not. fc, = Yl^—i A^ is 
the degree of the ith node and e is the overall coupling 
constant. In the present investigation we consider a ho- 
mogeneous delay r. The function f(x) defines a local 
nonlinear map, g(x) defines the nature of coupling be- 
tween the nodes. We consider phase synchronization as 
described in [l9|. As the network evolves, it splits into 
several synchronized clusters. In order to have a clear 
picture of SO and D behavior, we use fintra an( i finter 
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FIG. 1: Phase diagram of phase synchronization patterns in 
system (JXJ) for 1-d lattice with N = 50, < k >= 4. Gray-scale 
encoding represents values of (a) finter and (b) f intra Local 
dynamics is governed by logistic map f(x) = 4x(l — x) and 
coupling function g(x) — f(x). The figure is obtained by av- 
eraging over 20 random initial conditions. The regions, which 
are black in both graphs (a) and (b), correspond to states of 
no cluster formation. Both subfigures with gray shades corre- 
spond to clusters having both inter- and intra-couplings. The 
regions in (a), which are lighter as compared to the corre- 
sponding e and r values in (b), refer to dominant D phase 
synchronized clusters, and the reverse refer to dominant SO 
phase synchronized clusters. White regions in (a) and (b) re- 
fer to ideal D or ideal SO clusters respectively. The regions, 
which are dark gray in (a) and black in (b) or vice-versa, cor- 
respond to states where a much less clusters are formed, (c) 
and (d) are for scale-free, and (e) and (f ) are for bipartite net- 
works and demonstrate the same as (a) and (b) respectively. 



measures for intra- and inter-cluster couplings as follows 
[1; f intra = N intra /N c and finter = N mter /N c . where 
N c is total number of connections in the network, Ni ntra 
and Ni n t er are the numbers of intra- and inter-cluster 
couplings, respectively (20|. We evolve Eq.([T]) starting 
from random initial conditions, and study the dynam- 
ical behavior of nodes after an initial transient. First 
let us consider local dynamics being governed by the lo- 
gistic map f(x) = 4x(l — a;), and the coupling function 
g(x) = f(x). 

The undelaycd coupled maps on all model networks we 
have considered yield dominant D clusters in the range 
(0.16 < e < 0.25). For rest e values, coupled maps on 1-d 
lattice and small-world networks exhibit no phase syn- 
chronization, except for e > 0.74 having mixed clusters 
with very small values of finter and f lntr a (Figs. [J|a) 
and QIb)). In this e range scale- free and random net- 
works favor synchronization yielding better cluster for- 
mation than the corresponding regular and small-world 
networks (Figs.[]Jc) andQJd)), while bipartite networks 
lead to ideal SO synchronization for 0.45 < e < 0.85 and 
ideal D synchronization for higher e values 0] . Upon in- 
troducing a delay of r = 1 in Eq.([T]), after very small 
e values, for which there is no phase synchronization for 
the undelayed case (black color for the Figs.QJa), (b), (c) 
and (d)), we get SO clusters in the region 0.13 < e < 0.2 
as seen from the white regions in the Figs.rTJb) andQJd). 
For most of the e values in this region, coupled dynam- 



ics exhibits a periodic evolution with a period depending 
upon t. For a further increase in s, in the middle cou- 
pling range, the 1-d lattice, small-world, scale-free and 
random networks lead to an increase in D synchroniza- 
tion in 0.4 < e < 0.7, whereas for complete bipartite 
networks, ideal D synchronization is achieved for almost 
all e values in this range. For 0.85 < e < 1.0, the delayed 
case exhibits a very small (almost negligible) cluster for- 
mation compared to the undelayed case, hence indicat- 
ing a suppression of synchronization for all the networks 
except for bipartite networks forming ideal D clusters. 
For t = 2, lower e range coerces the formation of dom- 
inant D clusters, similar to the undelayed case. As e 
increases, 1-d lattice and small-world networks lead to 
mixed clusters, whereas scale-free and random networks 
lead to dominant D clusters. Bipartite networks emu- 
lates ideal D synchronization. With a further increase in 
t, at a lower e range, odd r leads to a similar behavior 
as for t = 1 and even r exhibits similar behavior as for 
t = and r = 2. For the intermediate e range there is a 
suppression in synchronization. Higher e values manifest 
no cluster formation as illustrated by the black regions in 
Fig. [1] for all networks except bipartite which form ideal 
D clusters for e > 0.4 for all r. 

Above description boils down to the following; there 
is a e region which demonstrates a change in the phe- 
nomenon of cluster formation with a change in r. The 
zero and even delays imply dominant D clusters, whereas 
odd delays imply ideal or dominant SO clusters. More- 
over, odd delays lead to SO clusters with a periodic evo- 
lution, whereas zero and even delays lead to D cluster 



with periodic, quasi-periodic or the chaotic evolution [21 
Note that the measure of phase-synchronization consid- 
ered here satisfies the metric properties, but does not in- 
clude anti-phase synchronization and consequently nodes 
being anti-phase synchronized would land up in different 
clusters. However, anti-phase or phase shift synchroniza- 
tion is not the only cause behind the separation of nodes 
in clusters [2lj . 

Though the nodes in various clusters display a rich 
dynamical evolution, a simple analysis for periodic syn- 
chronized state, for example bipartite networks in the 
lower e region, provides a basic understanding of differ- 
ent behaviors indicated by odd and even delays. In this 
e region, the coupling term having a delay part yields: 



/(*(* - r)) 



f(pi) if r = and even 
f{j>2) if T is odd, 



implying that the discrete time delay considered here 
introduces a difference on the evolution of the nodes 
(Eq.([T])) depending upon the parity of delay, and thus 
leading to a particular behavior for zero and even delays 
but a different behavior for odd delays. 

Furthermore, a change in r leads to a change in SO 
or D cluster pattern. A pattern refers to a particular 
phase synchronized state, containing information of all 
the pairs of the phase synchronized nodes distributed in 
the various clusters. A change in the pattern refers to 
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FIG. 2: (Color online) A typical behavior of coupled dy- 
namics illustrating D patterns observed with changing r. 
Squares represent clusters, diagonal dots represent freely- 
evolved nodes while off-diagonal dots imply that the two cor- 
responding nodes are coupled (i.e. Aij = f). In each case the 
node numbers are reorganized so that the nodes belonging to 
the same cluster are numbered consecutively. The example 
presents a scale-free network with N — 50 and e = 0.6. For 
r = 0, very few nodes are forming cluster. For r = 1, 3 and 5, 
nodes form dominant D clusters, whereas r = 2 and 4, yield 
very few nodes forming clusters of ideal D type. 



the state when members of a cluster get changed as an 
effect of delay. For some cases we observe ideal D or SO 
clusters. Ideal SO synchronization refers to a state when 
clusters do not have any connection outside the cluster, 
except one. The ideal D synchronization refers to the 
state when clusters do not have any connections within 
them, and all connections are outside. 

Next we focus to the e range where the delayed evolu- 
tion leads to ideal D clusters for bipartite networks, and 
dominant D clusters for other networks. In bipartite net- 
works, a division of nodes into ideal D clusters is unique. 
Whereas for other network structures, there can be var- 
ious possible ways in which one can distribute nodes to 
form ideal D (for average degree two) or dominant D (for 
larger average degree) clusters. Fig. [2] plots snapshots of 
clusters for different r by keeping all other parameters 
same. It indicates that with a change in t, both nodes 
forming clusters as well as size of clusters are changed. 
Note that the dynamical evolution here may be periodic, 
quasi-periodic or chaotic. In this region, for a particu- 
lar delay value, the clusters are almost stable with time 
evolution, with few nodes of the floating type Q. But 
a change of r has a drastic impact on cluster patterns, 
and may lead to entirely different sets of nodes form- 
ing clusters. Hence D patterns obtained in this range 
are dynamic with respect to a change in r. However, 
the phenomenon behind the pattern formation does not 
change, and the D mechanism is mainly responsible for 
the cluster formation. For this e range, a delayed evolu- 
tion on a bipartite network yields ideal D clusters for all 
t values we have investigated. 

Aforementioned can be explained further using the ex- 
ample of bipartite networks. A Lyapunov function anal- 
ysis can be carried out for delayed case in a very similar 
fashion as for r = described in [22j . and for a pair of 
synchronized nodes on a bipartite network can be written 



t-X 1 2 3 1.2.3 
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FIG. 3: (Color online) Three nodes schematic diagram illus- 
trating impact of delay. Arrows depict the direction of infor- 
mation flow as governed by Eq.JT]). The dashed lines show the 
flow of information from the (t — 2) th time step. For r = 0, 
evolution of all nodes (•) receive information from the second 
node (left panel), whereas in presence of delay, evolution of 
connected nodes at a particular time do not involve any com- 
mon term (right panel). For both panels, first and third nodes 
are connected with the second one leading to the construction 
of the smallest possible bipartite network. 



as: 



V ij (t + l) = [(l-e)(f(x i (t))-f(x j (t))) + 

2e N 2s N ' 2 

- J2 9(x j (t-r))--J29(xi(t-r))}' 



j=N/2+l 



For ideal D state, the synchronization between two nodes 
which arc not directly connected is independent of the de- 
lay terms as the coupling terms cancel out, and only de- 
pends on £. Hence, delay does not affect synchronization 
between the nodes which are not directly connected (22| , 
and only comprehends its presence for those which are 
directly connected. As a consequence, depending upon 
s and r, it may either enhance or destroy the synchrony 
between them. For instance, in the lower s range odd 
delays lead to an enhancement of coordination between 
connected nodes yielding a transition to SO clusters. 
Whereas in middle s range, delay destroys synchroniza- 
tion between connected nodes yielding D clusters state. 
As indicated by Fig. [3J for r = 0, the common term in 
the evolution equation for all the nodes may be the rea- 
son for global synchronization. Whereas, for r > 0, the 
network gets divided into two parts, one set of nodes has 
completely different terms in its evolution equations than 
those of the second set. An important inference of our 
results is that in the presence of delay, the dynamical 
evolution on bipartite network identifies the underlying 
network structure and gives rise to ideal D clusters for 
almost all the couplings for s > 0.4. Note that a previous 
result on delayed bipartite networks concludes that they 
would lead to worst synchronization @, but D clusters 
observed here very clearly reveal a very good synchroniz- 
ing power of the same. 

In order to demonstrate the robustness of the above 
phenomena, we also present results for coupled circle 
maps. In Eq.([T|), the local dynamics is defined by the 
circle map, f(x) = x + u> + (p/2ir)sin(2irx), with param- 
eter values taken in a chaotic regime. Fig. H plots the 
examples demonstrating the S-D transition, furthermore 
different r values are associated with a change cluster 
pattern as manifested by coupled logistic maps. 
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FIG. 4: (Color online) Phase synchronized patterns for cou- 
pled circle maps on scale-free networks with N — 50, < k >= 
2, g(x) = x and e = 0.24. 

Wc have studied effects of delay on phenomena of phase 
synchronized cluster formation in coupled map networks. 
Depending upon e values, a change in r may lead to a 
change in the phenomenon of cluster formation, with de- 
lays of the same parity being associated with the same 
phenomenon, or favors D clusters for delayed case over 
undclayed one which in extreme case of bipartite net- 
works demonstrate robustness of D mechanism against 
change in r. Furthermore, different r values may lead 
to an entirely new pattern of the cluster. For example, 
in middle e range, different r values lead to different dy- 
namical patterns of dominant D type. Whereas lower 
e values produce dynamical patterns of dominant D or 
dominant SO type. 



While an enhancement or suppression of complete syn- 
chronization as an introduction of delay was already well 
investigated in coupled maps models, mechanisms of de- 
layed unit to unit interaction were unknown. Delay may 
enhance the coordination among the connected nodes 
leading to an enhancement of synchronization identify- 
ing underlying connection topology, which had been the 
main theme of a few recent studies, but observation of 
a D mechanism behind the cluster formation in delayed 
coupled networks is a new insight suggesting that delay- 
induced synchronization may lead to a completely dif- 
ferent relation between functional clusters and topology, 
than relations observed for the undelayed evolution. Our 
study draws its significance in understanding synchro- 
nization in real world networks such as neural networks, 
where clusters are formed due to delayed interactions be- 
tween neurons [23| and may be of D type [24| . An analy- 
sis presented for bipartite and periodic cases help in dis- 
cerning a possible impact of r on the coupled evolution in 
such systems. Moreover, a change in patterns of neural 
activities has been found to be related with brain disor- 
ders such as Alzheimer (25^ . Research in the dimension 
of delay-induced patterns might propagate a finer appre- 
hension of the origin and treatment of these diseases. At 
fundamental level, a study of phase shift synchronization 
[H, based on phase synchronization measure considered 
here, is an aspect to explore in future [26j . 
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